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I. INTRODUCTION 

Masses of quarks and leptons are fundamental param- 
eters of the standard model. They cannot be determined 
by the theory and must be measured. A precise measure- 
ment of the mass of the r lepton is important for testing 
lepton universality [l[ and for calculating branching frac- 
tions that depend on the r mass • Uncertainties in the 
T mass have important consequences on the accuracy of 
the calculated leptonic-decay rate of the r, proportional 

to m5 g. 

CPT invariance is a fundamental symmetry of any 
local field theory, including the standard model. Any 
evidence of CPT violation would be evidence of local 
Lorentz violation and a sign of physics beyond the stan- 
dard model d, H, 0]- The most common tests of 
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CPT invariance are measurements of the differences of 
the masses and hfetimes of particles and their antiparti- 
cles. The most precise test of CPT invariance is from the 
measured hmits of the mass difference of neutral kaons, 
|M^o - M^\/M^y(. < 8 X 10-1^ 8] at 90% confidence 
level (CL), where M^yQ is the average value of M^o and 

At the T(4S') resonance, the cross section for e+e 
T^T^ is 0.919± 0.003 nb Q, resulting in a very large data 
sample, comparable to the number of bb events produced. 
With this data sample we can perform a pseudomass- 
endpoint measurement, first used by the ARGUS Col- 
laboration [l^ and recently by the Belle Collabora- 
tion [ll| , to measure the mass of the r lepton. Unlike the 
production-threshold method used by the BES T] and 
KEDR [12] experiments, this pseudomass method has 
the advantage of measuring the mass of the r+ and 
separately, which allows us to test the CPT theorem by 
measuring their mass difference. This measurement was 
first performed by the OPAL Collaboration , and the 
current limit is |M^+ - M^-|/M^^g < 2.8 x 10"'' H at 
90 % CL: the Particle Data Group (PDG) average value 
of the r mass is M^yQ = 1776.84±0.17 MeV Tables 
U and |TT] summarize the most recent measurements of Mr 
and the measured upper limits of \Mr+ — M^-l/M'^yQ. 




Pseudomass (GeV) 



FIG. 1: Pseudomass distribution. The points are data, the 
solid area is the background estimated from MC, and the 
dashed vertical line represents the PDG average value of the 
T mass 01 . Note the sharp edge of the distribution at the r 
mass. 

leptons. We also assume M^^ = 0. Since the neutrino 
is undetected, we cannot measure the angle 9*, thus we 
define the pseudomass Mp by setting 9* = 0: 



TABLE I: Recent r mass measurements. 



Experiment 



Mr (MeV) 



1776.96 



+0.18+0.25 



BES fl] 
KEDR [12] 
Belle [llj 1776.61 ± 0.13 ± 0.35 



1776.8ll^:23± 0.15 



TABLE IL Measured upper limits of the t'^ and r 
difference at 90% CL. 

Experiment |A/^+ - M^- \IM\yq 



OPAL 
Belle 



[13] 
[11] 



< 3.0 X IQ- 

< 2.8 X IQ- 



Ml 



2{^s/2-El){El-P*)<Mr. (2) 



Figure [T] shows the pseudomass distribution after ap- 
plying all of the selection criteria (Section IIII[) and the 
sharp kinematic cutoff at Mp = Mr- The smearing of the 
endpoint and large tail in the distribution is caused by 
ISR/FSR and detector resolution. The r mass is mea- 
sured by determining the position of the endpoint. We 
choose to use the decay mode t~ — > t:~ t:^ i:~ v^- and its 
charge conjugate, since it has a relatively large branch- 
ing ratio, B{t- ^ n-n+TT-j^r) = (8.99 ±0.06)% [|, has 
a high signal purity, and has large statistics in the end- 
point region due to the large central value and width of 
the mass distribution of the Stt system. 



The pseudomass is defined in terms of the mass, energy, 
and momenta of the r decay products. For hadronic 
decays of the (r^ h^Vr and its charge conjugate), 
the r mass, A/^, is given by 

Mr = ^Ml + 2(^/i/2 - i?*)(i?* - P* cos9*), (1) 

where 9* is the angle between the hadronic system and 
the I'r and Mh, Eh, and Ph are the mass, energy and mag- 
nitude of the three-momentum of the hadronic system 
h, respectively. The * indicates quantities in the e+e^ 
center-of-mass (CM) frame. In the CM frame, the energy 
of the r is given by E* = y/s/2, where i/s = 10.58 GeV. 
This relation ignores initial state radiation (ISR) from the 
e"'"e~ beams and final state radiation (FSR) from the r 



II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the 
BaBAR detector at the PEP-II asymmetric-energy e+e~ 
storage rings operating at the SLAG National Accelera- 
tor Laboratory. We use 423 fb^^ of data collected at the 
T(4S') resonance corresponding to over 388 million t+t~ 
pairs. For the control samples of inclusive A'^ — s- tt+tt^, 
D+ ~> K-'K+7:+, D+ ^ 07r+, D+ ^ 07r+, and their 
charge conjugates used for systematic studies, we use 
about 100, 100, 423 and 423 fb"^ of data, respectively. 
The background Monte Carlo (MC) samples used for 
this analysis comprise of generic e~^e~ T{AS)_—^ BB 
events simulated with the EvtGen generator [lj|, and 
e+e~ 99 (q = u, d, s, c) continuum events simu- 
lated with the Jetset7.4 generator [l^. For simulation 
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of r-pair events we use the MC generators KK2f flG*] and 
Tauola [13, and use PHOTOS Q^] to incorporate FSR. 
For the extraction of the t mass, we generate signal sam- 
ples with three different r masses {Mr = 1774, 1777, and 
1780 MeV), each comparable in event totals to the data 
sample. The BaBAR detector and its response to parti- 
cle interactions are modeled using the GEANT4 simulation 
package [Toj . 

The BaBAR detector is described in detail else- 
where |20|]. The momenta of the charged particles are 
measured with a combination of a five-layer silicon ver- 
tex tracker (SVT) and a 40-layer drift chamber (DCH) 
in a 1.5 T solenoidal magnetic field. A detector of inter- 
nally reflected Cherenkov radiation (DIRC) is used for 
charged particle identification. Kaons and protons are 
identified with likelihood ratios calculated from dE/dx 
measurements in the SVT and DCH, and from the ob- 
served pattern of Cherenkov light in the DIRC. A finely 
segmented CsI(Tl) electromagnetic calorimeter is used to 
detect and measure photons and neutral hadrons, and to 
identify electrons. The instrumented flux return contains 
resistive plate chambers and limited streamer tubes [2lj 
to identify muons and long-lived neutral hadrons. 

The most critical aspect of this analysis is the recon- 
struction of the charged particle momenta. Tracks are 
selected using the information collected by the SVT and 
DCH using a track finding algorithm: they are then refit 
usin^_a Kalman filter method to refine the track parame- 
ters 22]. This algorithm corrects for the energy loss and 
multiple scattering of the charged particles interacting 
with the detector material and for any inhomogeneities 
of the magnetic field according to a detailed model of 
the tracking environment. Since the energy loss depends 
on particle velocity, the Kalman filter is performed sep- 
arately for five mass hypotheses: electron, muon, pion, 
kaon, and proton. The main components of the detec- 
tor to be modeled for charged particle tracks originating 
from the vicinity of the interaction point are the 1.4 mm 
thick beryllium-beam pipe and 1.5 mm of cooling water 
at a radius of 2.5 cm, five layers of 300 fim thick silicon at 
radii of 3.3 cm to 15 cm, a 2 mm thick carbon-fiber tube 
at 22 cm that is used to support the SVT, and a 1 mm 
thick beryllium tube at 24 cm that makes up the inner 
wall of the DCH. Detailed knowledge of the material in 
the tracking volume and the magnetic field is crucial to 
accurate momentum reconstruction '^sj. This informa- 
tion is based on detailed information from engineering 
drawings and measurements taken both before and after 
the commissioning of the detector. 



III. ANALYSIS METHOD 

For our event selection, we require exactly four tracks 
in the event, none of which is identified as a charged kaon 
or proton. We veto events with K'^ tt+tt^ candidates 
with an invariant mass within ±25 MeV of the nominal 
mass [1] and tt" 77 candidates constructed with 



photons with CM energy greater than 30 MeV and an 
invariant mass in the range 100 MeV < M^j < 160 MeV. 
We require the total charge of the event to be zero. We 
divide the event into two hemispheres defined by the 
plane perpendicular to the event-thrust axis in the CM 
frame, which is calculated using all tracks and photon 
candidates. One hemisphere of the event, the tag side, 
must have a single track identified as either an electron 
or muon, and in the opposite hemisphere, the signal side, 
we require three charged tracks, none identified as a lep- 
ton. In addition to the tt^ veto, we require the number 
of photons with CM energy greater than 50 MeV on the 
signal side to be less than 5 and the total photon energy 
on the signal side to be less than 300 MeV to further 
reduce sources of background with one or more neutral 
pions. 

To reduce background events from two-photon pro- 
cesses, we apply six additional selection criteria. We re- 
quire the total reconstructed energy of the event to be 
within the range 2.5 GeV < E;^^ < 9.0 GeV and the 
thrust magnitude to be greater than 0.85, where these 
quantities are calculated with all tracks and neutrals with 
CM energy greater than 50 MeV. We require the tag 
lepton to have an energy less than 4.8 GeV, the en- 
ergy of the three pion system on the signal side to be 
1.0 GeV < E^^ < 5.2 GeV, and the reconstructed mass 
of the 37r system to be greater than 0.5 GeV. We also 
require the polar angle of the missing momentum to be 
in the range -0.95 < cos 9*^^^^ < 0.92. 

We define our fit region to be 1.68 < Mp < 1.86 
GeV. After all requirements, our signal efficiency is 2.0% 
and the purity of our sample is 96%. Our largest back- 
ground is —)■ 2Tr^ TT^ TT^ Ut , whcrc the ir'^ is not re- 
constructed. The total number of events in the data is 
341,614, 340,243, 352,609, and 329,248 for the r+, r", e 
tag, and ^ tag, respectively. 

We use three MC samples with different t masses 
(1774, 1777, and 1780 MeV) to empirically determine 
the relation between the pseudomass endpoint and the t 
mass, accounting for the smearing due to resolution and 
ISR/FSR effects. 

To determine the endpoint from the pseudomass distri- 
bution, we perform an unbinned-maximum-likelihood fit 
to the data using an empirical function ^11,] of the form 

F{x) = {p3 + pix) tan"^ ( — ) + + P&x, (3) 

where x is the pseudomass, and the pi are free param- 
eters of the fit. Only the position of the endpoint, pi, 
is important in determining the r mass, as the shape of 
the distribution does not affect the edge position since 
the correlation between pi and the other parameters is 
small. 

Figure [2] shows the pseudomass distributions from the 
three MC samples, with the shift in the edge clearly visi- 
ble. We fit each one of the MC distributions, and Figure 
|3]shows the fit results for pi versus the generated r mass. 
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FIG. 2: Pseudomass endpoint distributions from the three 
MC samples. The open squares, dots, and open triangles are 
for generated r mass values of 1774, 1777, and 1780 MeV, 
respectively. 

In the absence of ISR/FSR effects and with perfect de- 
tector resolution we would expect the relation between 
the pi fit result and generated t mass to be linear with a 
slope of unity and y-intercept = 0. With the inclusion of 
the ISR/FSR effects and detector resolution, we expect 
the relationship to still be linear with a slope of unity but 
to have a non-zero offset. We fit the results to a linear 
function, {pi - 1777 MeV) = ai(Mg - 1777 MeV) -I- oq, 
where Mg is the generated r mass, and uq and oi are 
free parameters of the fit. The results of the straight-line 
fit are ai = 0.96 ± 0.02 and ao = 1.49 ± 0.05 MeV. We 
use the results from the straight-line fit to determine the 
value of the r mass from the endpoint fit of the data. 

To determine the mass difference, we split our data 
sample into two sets based on the total charge of the 
three pion signal tracks. We use the combined fit results 
for oi and ao to determine the mass of t+ and t~ . As a 
cross check, we split our MC in the same way and repeat 
the procedure described above for each sample. We find 
the individual fit results for ai and ag to be within one 
sigma of the combined fit results. 



IV. TRACK MOMENTUM RECONSTRUCTION 

A previous analysis [2^ of BaBAR data has revealed 
that the track reconstruction procedure leads to system- 
atic underestimation of the individual track momentum. 
This effect is not observed for MC simulation. There 
are two potential sources of bias in the track momentum 
measurement: errors in the detector model which could 
lead to a momentum-dependent bias, and incorrect mod- 
eling of the magnetic field strength in the tracking vol- 
ume, which leads to a bias independent of momentum. 
We use a — > tt^tt control sample to investigate this 
bias and determine a correction. The K*^ daughter pi- 
ous have a momentum distribution similar to the pious 




2 ^ J 

-3-2-10123 

(M - 1777) (MeV) 



FIG. 3: The fitted value of pi as a function of Mg, the value 
of Mt in the simulation. These fit results are used in the 
determination of Mr from the endpoint fit to the data. 



in our signal sample, and the long fiight length of the 
Kg is ideal for studying the energy loss correction used 
by the reconstruction algorithm. The K'^ candidates are 
reconstructed from two oppositely charged tracks that 
have an invariant mass within 25 MeV of the nominal 
mass value [1]. This sample comprises 2.96 x 10^ K'^ 
candidates. We determine the mass by performing 
a maximum-likelihood fit to the data using a function 
which is a sum of two Gaussian distribution functions 
with a common mean and different widths, and a sec- 
ond order polynomial to describe the background. The 
background is relatively fiat and does not affect the mea- 
surement of the ifg mass. We increase the amount of 
SVT material, the strength of the solenoid field, and 
the strength of the field due to the magnetization of 
the beam-line dipole magnets inside the detector model 
to correct the reconstructed track momenta. The in- 
creases of the material in the SVT and the strength of the 
solenoid field are chosen to improve the agreement of the 
reconstructed K° mass with the world average value [1] . 
These increases are larger than the estimated uncertain- 
ties. In the following we detail the procedure to derive 
these corrections. 



A. Energy Loss 

Track momenta are corrected for energy loss by the 
Kalman filter procedure. The amount of energy a parti- 
cle loses due to material interactions depends on the na- 
ture and amount of material traversed and the type and 
momentum of the particle. Thus, any error in the esti- 
mated energy loss will vary with the amount of material 
the track traverses and the laboratory (lab) momentum 
of the track. 

There is clear evidence that track momenta are un- 
derestimated for our nominal reconstruction procedure, 
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as shown in Figure H The sample, as a function of 
the decay- vertex radius, ranges in purity from 7% to 91%, 
with the lowest purity arising from the interaction region, 
when the candidates have very short flight distances. The 
larger the radial distance of the if^ decay vertex, the less 
material the charged pions traverse, decreasing the size of 
the energy-loss correction. The largest deviation is seen 
for those events where the K^^ vertex is closest to the in- 
teraction region. This dependence of the reconstructed 

mass on the amount of material traversed by the pions 
demonstrates that the energy-loss correction is underesti- 
mated. Figure m also shows the mass as a function of 
the lab momentum: the purity of the sample ranges 
from 14% to 83% with increasing momenta. We see that 
lower momenta i^g particles have masses further from 
the expected mass than high momenta ones, since 
the energy-loss corrections are greater for the lower mo- 
menta particles, because the -fCg decay-vertex radius and 

momentum are correlated. 
We study two possible corrections to the energy-loss 
underestimation: increasing the amount of SVT mate- 
rial by 20% and increasing the amount of material in 
the entire-tracking volume by 10% [2^. For each cor- 
rection, we increase the density of the corresponding de- 
tector material by the indicated amount and repeat the 
Kalman filter procedure again. Figures [Hand [5] show the 
resulting if" mass variations after these corrections. In 
the case where the entire-tracking-volume is increased, 
we observe that the K^^ mass variation with the decay- 
vertex radius is flat, but the mass is over-corrected 
at lower momenta. A smaller correction of the entire 
tracking material could be used to flatten the mass 
variation with the momentum, but then the mass 
variation with the decay-vertex radius would no longer 
be flat. Therefore, we do not use the increase of the en- 
tire tracking material in our correction. In the case where 
the SVT material is increased, we observe that the ifg 
mass variation with decay-vertex radius and momentum 
is substantially reduced and flat. This reduction and flat- 
tening of the dependence of the reconstructed mass 
on the vertex position and momentum is our motivation 
for applying this correction to our material model. The 
estimated uncertainty in the SVT material is about 4.5% 
as determined from detailed analyses of the composition 
of the SVT and its electronics:, the 20% increase signif- 
icantly exceeds this estimated uncertainty. We apply it 
as a simplified method to account for this and all other 
uncertainties in the energy loss estimation. The uncer- 
tainty in this simple correction accounts for the largest 
uncertainty in the r mass measurement. 



B. Magnetic Field 

After the SVT energy-loss correction, we find that the 
ifg mass is still underestimated. In order to further cor- 
rect the mass, we consider two possible sources of 
error: uncertainty in the 1.5 T solenoidal magnetic field 
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FIG. 4: Fitted ifg mass vs. decay-vertex radius (top) and 
i^s lab momentum (bottom). On the vertical axis, the PDG 
average value of the K'l mass 0| has been subtracted from 
the fitted value. The points show the normally reconstructed 
data events, the open circles show the data reconstructed with 
20% more SVT material, and the shaded region is the error on 
the nominal K% mass ! 8j . The dependence of the K% mass on 
the decay-vertex radius and momentum is due to the under- 
estimation of the energy loss by the reconstruction procedure. 



that runs parallel to the beam axis and the perturbation 
to this field due to the magnetization of the magnetic ma- 
terials comprising the beam-line dipolc magnets (BDM) 
due to the solenoid field. 

The BDM are permanent magnets, made of samarium- 
cobalt, the closest of which is located 20 cm away from 
the interaction region. The fringe fields from these mag- 
nets in the interaction region are small and have been 
well measured; however the magnetic field due to the 
magnetization of these magnets by the solenoid field is 
not well known. The permeability of the BDM mate- 
rial was not measured before the commissioning of the 
detector, and subsequently variations in the susceptibil- 
ity of ± 20% with respect to the average value (-1-0.14) 
have been found within the individual small blocks used 
to construct the BDM. The field in the tracking volume 
due to the magnetization was estimated from measure- 
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FIG. 5: Fitted Kg mass vs. decay-vertex radius (top) and lab 
momentum (bottom). On tlie vertical axis, the PDG average 
value of the Kg mass has been subtracted from the fit- 
ted value. The points show the normally reconstructed data 
events, the squares show the data reconstructed with 10% 
more material in the tracking volumCjand the shaded region 
is the error on the nominal if? mass 0. 



FIG. 6: Fitted K^ mass vs lab momentum. On the verti- 
cal axis, the PDG average value of the K^ mass has been 
subtracted from the fitted value. The points show the nor- 
mally reconstructed data, the open squares show the data 
reconstructed with the solenoid field increased by 0.02%, the 
crosses show the data reconstructed with magnetization field 
increased by 20%, and the shaded region is the error on the 
nominal Kq mass [al. 



ther correct the K'^ mass, we increase the solenoid field 
by 0.02%, and then refit the tracks. Figure [6] shows the 
effect of this increase. This increase is larger than the 
measured uncertainty in the solenoid field, but it further 
improves the agreement of our reconstructed mass 
with the world average value [8]. Table Iml shows that 
this increase, in conjunction with the increases of the 
SVT material and the BDM magnetization field, shifts 
the ifg mass so that it is consistent with the world aver- 
age value [1]. 



Momentum Reconstruction Correction 



ments made at two points near the BDM, using Hall and 
nuclear magnetic resonance probes, followed by finite el- 
ement calculations that depend on the permeability of 
the magnets. In 2002, the probes were moved and the 
field was re-measured at two new points. At one point, 
there was good agreement with expectation, but at the 
other point the overall value of the field strength was 
0.4% higher than expected. We increase the field due to 
the magnetization of the BDM by 20% to account for the 
variation of the measured permeablility of these magnets 
and the observed discrepancy between the measured and 
estimated fields, in order to improve the agreement of our 
reconstructed K'^ mass with the world average value [1] . 
Figure [S] shows the effect of the increase on the K'^ mass 
as a function of the momentum measured in the lab 
frame. 

The solenoid field was very accurately measured with 
an uncertainty of 0.2 mT prior to the installation of the 
BDM during the commissioning of the detector. To fur- 



We study the overall corrections described above for 
and decays. These corrections affect the recon- 
structed masses in different ways. Although the size of 
the correction varies depending on the decay kinematics, 
decay mode, and the mass of the particle being recon- 
structed, the masses of our test samples are consistent 
with the world averages after the corrections are applied. 
To determine the size of a correction, we increase the 
amount of SVT material, the field due to the BDM, and 
the solenoid field strength in three separate simulations. 
For each of these simulations we refit each pion track 
using the Kalman fit procedure described above, and re- 
calculate the reconstructed mass. The overall mass cor- 
rection is taken as the sum of the three individual mass 
shifts, and the corrected mass is determined by adding 
the correction to the mass determined with the normal 
reconstruction. Figure [7] shows the corrected K'^ mass as 
a function of the K'^ decay- vertex radius and momentum 
in the lab frame. This method improves the agreement 
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of our measured K'^ mass with the world average. Table 
mil shows the individual corrections as well as the overall 
correction for the mass. 

We also apply this method to the decay D'^ 
K~tt^t:^ and its charge conjugate. We perform a vertex 
fit to the three tracks and require the vertex probabil- 
ity to be greater than 0.1%. We also require the mass 
of the candidate to be in the range 1.84 GeV < AIo < 
1.90 GeV. To determine the mass of the D meson, we 
perform a maximum-likelihood fit to the Kittt mass dis- 
tribution using a function which is a sum of two Gaussian 
distribution functions with a common mean and differ- 
ent widths and a second order polynomial to describe 
the background. Table IIVI summarizes the result of the 
fits with the normal reconstruction and modified detec- 
tor model. We find that the measured mass using the 
normal reconstruction differs by —0.92 MeV relative to 
the world average value of 1869.62 ± 0.20 MeV Q; af- 
ter applying the correction, the difference is reduced to 
-1-0.15 MeV, in very good agreement within the current 
uncertainties. 

We apply this method to the events in the sample of 
T~ — !■ TT^TT+TT^i^^ and its charge conjugate and obtain a 
correction of -1-0.63 MeV for the r mass. Table fVl shows 
the individual shifts on the r mass. 



TABLE III: Shifts in the measured mean value of the Kg 
mass when each track reconstruction correction is applied sep- 
arately and comparison with the nominal value [1]. 



Fit 


M^o (MeV) 


Mass Shift ( MeV) 


Default Reconstruction 


497.323 




SVT Material +20.0% 


497.477 


4-0.154 


Solenoid Field +0.02% 


497.383 


+0.060 


BDM Field -H20.0% 


497.382 


+0.059 


Fully Corrected 


497.596 ± 0.006 




PDG Average 


497.614 ± 0.024 
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FIG. 7: Fitted mass vs. decay-vertex radius (top) and lab 
momentum (bottom). On the vertical axis, the PDG average 
value of the Kg mass @] has been subtracted from the fitted 
value. The points show the normally reconstructed data, the 
triangles show the data after the correction from the increased 
material and magnetic field strengths has been applied, and 
the shaded region is the error on the nominal Kg mass 



TABLE IV: Shifts in the measured mean value of the 
mass when each track reconstruction correction is applied sep- 
arately and comparison with the nominal value 



Fit 


Mo± (MeV) 


Mass Shift ( MeV) 


Default Reconstruction 


1868.70 




SVT Material +20.0% 


1869.17 


+0.47 


Solenoid Field +0.02% 


1869.00 


+0.30 


BDM Field +20.0% 


1869.00 


+0.30 


Fully Corrected 


1869.77 ± 0.04 




PDG Average 


1869.62 ± 0.20 





D. Charge Asymmetry 



The tt"*" and tt tracking efficiencies differ because of 
different cross sections for interactions of low-momentum 



TABLE V: Observed shifts for Mr in the data due to each 
correction applied to the reconstructed track momenta sepa- 



Detector Parameter 


Mr Shift ( MeV) 


SVT Material +20.0% 


+0.31 


Solenoid Field +0.02% 


+0.11 


BDM Field +20.0% 


+0.21 


Correction 


+0.63 



7r+ and tt^ with the detector material [8|. This could 
cause differences between the reconstruction efficiencies 
for and — > 'k~ tt'^ tt'^ and t~ — s- tt^ tt~ . A 
difference in the reconstruction efficiency for the t+ and 
might introduce a dependence of the reconstructed r 
mass on the r momentum and thus might result in an 
artificial mass difference. 

To estimate any charge asymmetry in the track recon- 
struction procedure, we measure the mass difference in 
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several control samples: K~TT'^Tr~^, _D+ — > 4>Tr~^, 

Df — > ^TT"*", and their charge conjugates. The momen- 
tum spectra of the daughter pions in the r signal sam- 
ple are similar to the spectra in the three control sam- 
ples. The selection criteria for the Z?"*" K~Tr'^TT~^ and 
charge conjugate modes are described in Section IIV CI 
The (f) candidates are reconstructed from two oppositely 
charged kaons, and the reconstructed mass of the 4> can- 
didate is required to be within ±12 MeV of the nomi- 
nal value Q. To reconstruct a D or Dg candidate, the 
two kaon tracks from the (j) candidate are combined with 
a pion track, and the three tracks are required to have 
a vertex probability greater than 0.1%. The D and Ds 
candidates are required to have a CM momentum greater 
than 2.4 GeV to further reduce backgrounds. The D and 
Ds candidates are required to have an invariant mass 
within the range 1.85 GeV < Md < 1.90 GeV and 
1.95 GeV < Mds < 1-99 GeV. For the three samples 
respectively, the total numbers of events arc 4.5 x 10^, 
1.7 X 10^ and 2.2 x 10^, and the purities are 33%, 
90%, and 87%. To determine the masses, we perform a 
maximum-likelihood fit to each three-particle invariant- 
mass distribution, again using a sum of two Gaussian 
distribution functions with a common mean and differ- 
ent widths and a second order polynomial background 
function. Table IVII shows the observed mass difference, 
AM = Mx+ — Mx- , for each of the three decay modes, 
where X is the particle whose mass is measured. The 
results are consistent with zero difference. Thus, we do 
not make any correction and use these results to deter- 
mine the systematic uncertainty in — M^- due to 
the possible residual uncertainty in tracking. As a cross 
check, we perform the study on a sample of D"*" — *■ (/'tt'^, 
— > (pTT'^ , and charge conjugates where we do not con- 
strain the momentum of the D and Ds ■ We find the mass 
difference of these samples is consistent with the results 
using the samples that have a D and Ds momentum con- 
straint. 



Mr = 1776.68 ± 0.12 MeV, where the error is statisti- 
cal only. 

Figure [5] shows the resulting pseudomass distribution 
from subtracting the t~ distribution from the t'^ distri- 
bution. We measure pi(r+) to be 1777.29 ±0.16 MeV 
andpi(r-) = 1777.88 ±0.17 MeV. Applying the above 
procedure, we find Af^+ — 1776.38 ± 0.16 MeV and 
Mr- = 1776.99 ± 0.17 MeV, where the errors are statis- 
tical only. Thus, Mr+-Mr- = -0.61±0.23(stai) MeV. 
Figure [TO] shows the pseudomass threshold region, where 
the T~ distribution is clearly shifted to a higher mass 
relative to the r"*" distribution. 
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FIG. 8: Combined r and pseudomass endpoint distri- 
bution. The points show the data, the curve is the fit to 
the data, and the solid area is the background. The inset is 
an enlargement of the boxed region around the edge position 
showing the fit quality where p\ is most sensitive. 



TABLE VL AM for the and Df meson control samples 
used to study the possible charge asymmetry. 



Sample 


Mass Difference ( MeV) 


D+ -» K~-K + Tl + 


-0.04 ± 0.03 


D+ -> (j)TT + 


±0.06 ± 0.04 


D+ ^ ,^7r+ 


±0.10 ±0.05 


V. 


RESULTS 



Figure [8] shows the pseudomass distribution of the 
combined t+ and t~ samples compared to the fitted 
distribution. The fitted value of the endpoint position 
is pi = 1777.58 ± 0.12 MeV. Using the MC resuhs 
for oq and Oi and applying the reconstruction proce- 
dure corrections described in Section IIV C[ we obtain 
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FIG. 9: Resulting pseudomass distribution from subtracting 
the T~ distribution from the distribution. 
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FIG. 10: Pseudomass distributions for the r+ and t~ in the 
region around the pseudomass threshold region. The open 
circles and solid points show the r+ and r~ distributions, 
respectively. The curves show the results of the fits to the 
data. 



VI. SYSTEMATIC STUDIES 

Table IVlIl summarizes the estimated systematic uncer- 
tainties in Mt- 

The largest source of error in the r mass measurement 
arises from the momentum-reconstruction uncertainty. 
The increases in the SVT material and magnetic-field 
strengths are applied to obtain a better agreement of the 
reconstructed K'^ mass with the nominal if^ mass Q , but 
the actual cause of the discrepancy is still unknown. The 
effect of the induced magnetization of the BDM on the 
magnetic field in the tracking volume has never been mea- 
sured, and the discrepancy between the actual field and 
modeled field is unknown. Although there is no evidence 
of any mis-modeling of the solenoid field, we increase the 
field by 0.02%, which is larger than the measured uncer- 
tainty of the field, 0.013% (0.2 mT). The simulation of 
the SVT material is believed to be accurate to within 
4.5%, but the increase we use is substantially larger. To 
account for the uncertainty of the momentum reconstruc- 
tion we add the mass shifts originating from these cor- 
rections in quadrature. This results in the dominating 
systematic uncertainty of ±0.39 MeV. The systematic 
uncertainties are summarized in Table |Vl 



TABLE VII: Systematic uncertainties in Mr- 



Source 


Uncertainty ( MeV) 


Momentum Reconstruction 


0.39 


CM Energy 


0.09 


MC Modeling 


0.05 


MC Statistics 


0.05 


Fit Range 


0.05 


Parameterization 


0.03 


Total 


0.41 



Another important source of systematic error comes 
from the uncertainty in the absolute scale of the e'^e~ 
CM energy. From the error propagation of Equation [21 
we find 

Near the endpoint of the pseudomass distribution, « 
y/s/2, and M* « M^, so that 6{Mr) w 0.17<5(Vs/2). 

The e+e" CM energy calibration has been seen to drift 
over time due to changing beam conditions. Over a two 
year period of data taking, the calibration had drifted by 
—2.6 MeV. We exploit the fact that the T{AS) resonance 
decays exclusively to bb pairs and calibrate Vs/2 based 
on the measured invariant mass of fully reconstructed B 
meson decays using the equation 

Mb = ^{V~s/2y~P*^ (5) 

where Mb and Pb are the mass and reconstructed mo- 
mentum of the B meson. We reconstruct a dozen 
hadronic B decay modes and divide the data into sub- 
samples of 2500 candidates each. We then perform a 
maximum likelihood fit to the reconstructed mass distri- 
bution for each subsample to extract the central value of 
Mb and then adjust to obtain the world average 

B meson mass [8[ . We apply this correction to the value 
of \/s/2 for all data taken during the time period corre- 
sponding to each subsample. The statistical uncertainty 
of this correction is negligible, so the only uncertainty 
in -y/s/2 is due to the error in the PDG average value 
of the B meson mass (0.5 MeV) 8]. This uncertainty in 
\/s/2 corresponds to a systematic uncertainty in M^ of 
0.09 MeV. 

Since we have a limited number of MC events, there 
are statistical errors associated with the straight-line fit 
parameters oq and ai (Figure [3]). These errors introduce 
a systematic error in Mr of ±0.05 MeV. 

We also consider alternatives for the pseudomass fit 
parameterization (Equation [3]), by fitting with two other 
functions [ll|: 



Fi{Mp) = {p3+p4Mp)— ^/ = +p5+P6Mp 

(6) 

and 

F2(Mp) = (p3 +P4Mp) ~^M,^p, +P5 +PeMp. (7) 



We repeat the fitting procedure with Fi and F2 and ob- 
tain shifts in Mr of -0.02 MeV and +0.02 MeV, re- 
spectively. We add the shifts in quadrature and find 
±0.03 MeV as the systematic uncertainty. 
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We also investigate the choice of fit range. We applied 
the procedure discussed in Section |TTT] using toy MC sam- 
ples, refitting each sample with various fit ranges. We 
take the largest shift, 0.05 MeV, as the systematic un- 
certainty. 

We study the effect of the MC modeling of the three- 
pion mass distribution in tau decays. We find that 
the peak of the distribution in MC is about 300 MeV 
lower than that in the data, while the widths of the 
distributions arc similar. The MC modeling for the 
TT^TT+TT^i'^ and its charge conjugates is based on 
16 form factors [2j| determined from low statistics data 
from the LEP and CLEO experiments: measuring the 
form factors is a very challenging task that has not yet 
been performed on the high statistics data collected by 
BaBAR. Although there is this discrepancy, we find that 
the pseudomass distribution in MC is similar to that in 
data. To test for possible effects in the endpoint of the 
pseudomass distribution due to the modeling of the Stt in- 
variant mass, we generate four toy MC samples, varying 
the mean and width of the Stt mass by ±300 MeV. We 
find that the shifts in the pseudomass endpoint are con- 
sistent with zero, but we conservatively take the average 
of these shifts, 0.05 MeV, as the systematic uncertainty 
due to the MC modehng. 

We also investigate the choice of background estima- 
tion and pion misidentification and find the effects on the 
fit result are negligible. We also find the error due to the 
uncertainty in the boost of the CM frame and the uncer- 
tainty in the MC modeling of the track resolution to be 
negfigible. 

We have assumed that the neutrino mass is zero even 
though the PDG limit for the direct measurement is 
M^^ < 18.2 MeV [Sj. Neutrino experiments [H] have 
measured differences in the mass squared between the 
three neutrinos to be much less than 1 eV^ [l^. Di- 
rect measurements of M^^ < 2 eV [1] thus suggest that 
the mass of the r neutrino is 0{< leV). We perform 
MC studies on the effect of the neutrino mass on the r 
mass determination and find that a 1 MeV neutrino mass 
would bias our result by -0.02 MeV. 

All of the systematic effects listed above cancel in the 
and T~ mass- difference measurement. An additional 
systematic arises from the possible charge asymmetry 
discussed in Section flVDI To study this effect, we mea- 
sure the mass differences for charged D and Ds mesons, 
which are presented in Table IVII We take a weighted 
average of the absolute values of the mass differences, 
0.06 MeV, as the resulting systematic uncertainty. As a 
cross check of the t sample, we studied the mass differ- 
ence Mt+ — Mt- separately for the e and /i tags, before 
and after the 20% increase of the SVT material, and find 
consistent results. 



VII. CONCLUSIONS 

In summary, we have measured the mass of the tau lep- 
ton to be 1776.68 ± 0.12(siai) ± GAl{syst) MeV, where 
the main source of uncertainty originates from the uncer- 
tainty in the reconstruction of charged particle momenta. 
This result is in agreement with the world average Q. 

We measure the mass difference of the t+ and t~ 
to be -0.61 ± Q.2i{stat) ± QM^syst) MeV, or (M^+ - 
Mr-)/M\y(j = (-3.4±1.3(sta<)±0.3(syst))xl0-'^. We 
use our result to calculate an upper limit on the mass dif- 
ference, \M^+-M^-^ M^Iavg < 5.5 X 10-4 at 90% CL. We 
find our measurement is consistent with the previously 
published results made by the Belle Collaboration. We 
perform parameterized MC studies to determine the sig- 
nificance of our result of the mass difference. We generate 
4500 samples each for the t+ and with the masses of 
each sample set to the value extracted from the combined 
data sample, 1776.68 MeV. The samples are generated 
with the same number of events as the number of events 
in the data. We fit each sample and calculate the mass 
difference between the t+ and samples. We also re- 
peat the procedure using an alternative parameterization 
(Equation [5]) , and determine that the two parameteriza- 
tions give consistent results. We find, assuming no CPT 
violation, that there is a 1.2% chance of obtaining a result 
as different from zero as our result. 
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